Dynamical restoration of symmetry breaking in realistic optical lattices 
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We show that properties of bosons loaded into the p-band of a rectangular optical lattice are 
substantially affected by the anharmonicity of the lattice potential. Even for very deep lattices, the 
optical confinement treated as a harmonic potential is highly oversimplified and in the low tunneling 
limit it leads incorrectly to a degeneracy of the ground state. We show that in this limit the degener- 
acy is lifted by the anharmonicity. For finite range of larger tunnelings the degeneracy is dynamically 
restored and the system gains an additional symmetry. This symmetry is, however, spontaneously 
broken in the thermodynamic limit and the state with non vanishing staggered angular momentum 
becomes the true ground state of the system. 



Confining ultra-cold gases in the optical lattices has be- 
come a standard experimental tool for studying strongly 
correlated many-body states [l|-l3[. In the last decade it 
became possible to arrange experiments in such a way 
that they can mimic different 'toy-models' known from 
condensed-matter physics. In this way, optical lattices 
are truly becoming dedicated quantum simulators for 
the study of different many-body configurations |^]. In 
this context, there is a recent interest to study ultra- 
cold atoms trapped in higher bands of optical lattices. 
Such states can be prepared by loading ultracold atoms 
to the higher bands_|5|-|8| or via interaction-induced inter- 
orbital transfer |9l4l3l|. Atoms in such excited bands can 
give rise to a zoo of exotic phases like, exotic superflui dity 
with broken time- reversal symmetry (TRS) |6l49l ll4l4l8l |. 
dynamical topological insulators [12'|, flat-band crystal- 
lization [19| etc. 

Recently, a speciflc model for p-orbital bosons in an 
effective one-dimensional optical lattice was introduced 
[20| . The authors assumed a local isotropy of lattice sites 
and high anisotropy of hopping amplitudes. In the limit 
of vanishing tunneling, they found that the system shows 
anti-ferro-orbital phase TRS breaking for the filling fac- 
tor V — 3/2. That is due to the fact that the local on-site 
Hamiltonian commutes with the angular momentum op- 
erator within harmonic approximation (HA) for the lat- 
tice sites. But recently it was shown in the context of 
Einstein-de Haas effect .ISj that anharmonicity is crucial 
when one study p-band physics. Taking into account the 
realistic optical potential results in an on-site Hamilto- 
nian which does not commute with the angular momen- 
tum. Motivated by this observation we have studied the 
system as described in |20| by taking into account the full 
complexity of the problem. We calculate all parameters 
of the Hamiltonian using the adequate single particle ba- 
sis of Wannier functions. We flnd that the degeneracy of 



the ground state which was destroyed by the anharmonic- 
ity in the low tunneling limit is surprisingly restored in 
the flnite range of larger tunnelings. Consequently, the 
possibility of TRS breaking states is investigated. 

The starting point of our analysis is the most general 
Hamiltonian describing spineless bosons conflned in the 
optical lattice potential Viatt('') interacting via contact 
interactions. The Hamiltonian has a form: 



ii = /dr*t(r) 



|;;V» + K.„,.) 



^(r) 
/dr^t(r)v^t(^)^(r)^(r)_ (1) 



The field operator ^(r) annihilates a boson at point 
r. Parameter g describes strength of the contact inter- 
actions and it is proportional to the s-wave scattering 
length. We expand the held operator ^{r) in single- 
particle basis of maximally localized Wannier states 
0f(r): 



*w = EE«««<^"(^)- 



(2) 



The operator da{i) annihilates a single boson at site i 
occupying Bloch band a. of the periodic potential. In the 
case of a rectangular two-dimensional optical lattice, i.e. 
when 

T4att(r) = 14 s\n^{2TTx/a.^) + Vy sin^ {2TTy / ay) + — ^z^ 

(3) 
the single-particle Hamiltonian is separable and Wannier 
states are products of one-dimensional Wannier functions 
4>t{r) = Wr/(x)Wf/(2/)Zo(z), where i = {i,,iy) and 
a — {ax, ay). We assume that the dynamics in the z 
direction is completely frozen and all bosons occupy the 
lowest state of the harmonic confinement Zo{z). 



The standard Bose-Hubbard Hamiltonian describing 
the lowest band dynamics is obtained by restricting the 
decomposition ^ to the lowest band only, i.e. a 



ing Hubbard-like Hamiltonian describing dynamics in ID 
chain: 



= 0. Such a Hamiltonian was extensively studied and "^ ^ X! ^^^^ ^ X! [^x^li'^)^^^) + tvO-li'^)^yU)] ■ (4a) 



many of its properties were discussed and verified exper- 
imentally [2l|. In this paper we are interested in many- 
body properties of bosons occupying higher orbitals of 
optical lattice. Therefore, we assume that the system is 
prepared in such a away that all particles occupy only the 
first excited Bloch band (p-band) of the optical lattice. It 
is known that in the case of non-interacting bosons in the 
square lattice (i.e. when Vx = Vy and a^ — ay) p-orbital 
states are degenerate and therefore many interesting pos- 
sibilities become open. Here we explore another interest- 
ing case of highly non-symmetric lattice with Vy ^ Vx- 
We choose Vy large enough to suppress the tunneling pro- 
cesses in y direction. At the same time we change the 
ratio ax I ay between lattice constants to preserve single- 
particle degeneracy between considered orbitals. In this 
way we end up with the system which can be effectively 
treated as a one-dimensional chain. At each site par- 
ticles can occupy two degenerate single-particle orbitals 
Pa; and Py. They can also tunnel to neighboring sites 
in X direction: yx (Py) orbital tunnels with amplitude 
tx < [ty > 0) and \tx\ > \ty\. The difference in sign 
and magnitude of the tunneling amplitudes is a direct 
consequence of the fact that Py {px) orbital tunnels with 
respect to the ground (excited) Bloch band of the opti- 
cal lattice in x direction. An oversimplified case of such 
a scenario was discussed originally in [20| . It was shown 
that for small tunneling the ground state of the system is 
doubly degenerate and in the thermodynamic limit the 
system undergoes spontaneous symmetry breaking. This 
result was obtained under an important assumption that 
each lattice site may be treated as a harmonic one. In a 
consequence, the single-particle part of the Hamiltonian 
becomes rotationally invariant because in this particu- 
lar approximation contact interactions between bosons 
preserve the local rotational symmetry. This is no longer 
true when one takes into account the anharmonicity of lo- 
cal sites. Here we reexamine properties of such a system 
by taking exactly the assumed geometry of the lattice. 
We find that properties of the system crucially depend 
on the anharmonicity and the resulting anisotropy of lat- 
tice sites. 

From now on we assume that bosons can occupy only 
the p-band states of the one-dimensional optical chain 
and we restrict the decomposition ^ to ct — (1,0) 
and ct — (0,1). Just for simplicity we introduce a new 
notation for corresponding annihilation operators ax{i), 
ay{i). These operators annihilate boson in the single- 
particle Wannier states (jjfir) = wf ^(x)>vf ^(2/)Zo(z) 
and (l)^{r) = W,} {x)Wl {y)Zo{z) respectively. We also 
introduce density operators for each orbital fiaij) — 
al^{j)aa{j) (a = x,y). In this way we obtain the foUow- 



The local, on-site Hamiltonian H{j) has the form: 



H{J) - E 
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All C/'s denote contact interactions between appropriate 
orbitals. Ex and Ey are single particle energies which 
in general are different. It is obvious that the Hamilto- 
nian commutes with the total number of particles oper- 
ator N = Nx + Ny, where Na — J2i "-q(*)j the property 
not enjoyed by Nx and Ny separately, due to the last 
two terms in the local Hamiltonian that transfer pairs of 
bosons between different orbitals. Thus the Hamiltonian 
has additional global Z2 symmetry related to the parity 
of the operator Ny (choosing Nx leads to the same con- 
clusions) and it commutes with the symmetry operator 
S — exp{iTrNy). To find the ground state of the sys- 
tem one can find the lowest energy states in two eigen- 
subspaces of S independently. Let us call these states 
\Geven) and \Godd) with corresponding eigenenergies Eeven 
and Eodd (subscripts even/odd correspond to even/odd 
number of bosons in orbital y). Finally one can choose 
the state with lower energy as the global ground state 
(GS) of the system. In principle it may happen that both 
ground states have the same energy. In such a case any 
superposition cos(0)|Gei,en) + sin(0)e"^|Godd) is a ground 
state of the system. As explained later, in the thermody- 
namic limit this U{1) x C/(l) symmetry is spontaneously 
broken to Ising-like Z2 symmetry and only one of two 
macroscopic states can be realized. 

Now, let us discuss the role of anharmonicity of the 
lattice potential. To make the discussion as clear as pos- 
sible we revisit the situation discussed in [20| and we 
take lattice depths Vx/Er^x — 6, Vy/Ea^y = 24 and fill- 
ing V = 3/2, where E^.a = 2'K^h /ma^. The single- 
particle energies Ex and Ey^ even for different lattice 
depths, can be equalized by changing the lattice con- 
stants ax and ay. For chosen Vx,Vy and in HA one can 
show that ax/ ay = 1/a/2 leads to equal single-particle en- 
ergies. Nevertheless, the ratio ax /ay calculated directly 
in the basis of Wannier functions differs from that value 
(see Fig. [TJa)) and for the example studied it is approx- 
imately equal to 0.42. This is the first place where the 
anharmonicity of the lattice potential is important. 

Importantly the anharmonicity dramatically changes 
also the contact interactions. Since the wave-functions of 
p-orbitals are just products of one-dimensional Wannier 
functions the ratios Uyy/Uxx and Uxy/Uxx do not depend 
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FIG. 1: (color on-line) Influence of the anharmonicity on the 
parameters of the Hamiltonian. (a) The ratio ax/ay in the 
HA (dotted blue line) and for Wannier basis (solid red line) 
equalizing single-particle energies E^ = Ey. (b) The ratio 
Uxx/Uyy (thick red line) and 3Uxy/Uyy (thin black line) cal- 
culated for exact Wannier functions. In both plots we take 
Vy/En^y = 24. The dotted vertical line indicates the case 
Vx/ Er^x ~ 6 studied in text. In HA both quantities are equal 
to one. 



on the lattice constants a^ and ay. Therefore they are 
functions of dimensionless lattice depths Vx/Efi^^x, and 
Vy/Ej^^y only. Moreover, in the HA they are equal to 
1 and 1/3 respectively, independently on lattice depths. 
This simple observation indicates that HA may be valid 
only in a very deep lattices regime. It is straightforward 
to show, that in the HA the Hamiltonian (|3]) reduces to 
the Hamiltonian which preserves angular momentum in 
each lattice site independently. The precise values of ra- 
tios Uyy/Uxx and Uxy/Uxx calculated from appropriate 
Wannier wave functions are presented in Fig. lb. We 
see that whenever dimensionless lattice depths are dif- 
ferent one gets Uxx 7^ Uyy. In the case studied in [20| 
the difference is larger than 10%. Moreover, even for 
equal lattice depths the contact interaction Uxy is never 
equal to Uxx/i- It means that harmonic limit can not 
be reached in any realistic optical lattice and rotational 
invariance of local lattice sites does not hold. 

The violation of the local rotational invariance is 
clearly visible when we change the single-particle basis 
from the Cartesian to the angular one. By introducing 
angular- momentum-like annihilation operators a±{j) = 
[ax{j) ± iay{j)\ /\/2 the local part of the Hamiltonian @ 
can be written in the form 



m-"^ 



n{j)[n{j)~'^^-\Ll{]) 



+ A 



{h{j)-l)[L+{j)+L^{j) 
\lI{j)-2.(l+{j)-L^{j)] -n(j) 



(5) 



where U = {Uxx + Uyy)/2, 5 = {Uxx - Uyy)/2, and A = 
Uxy — U/3. For convenience we introduce local density 
operator h{j) — a_^{j)a+{j) + a_{j)a-{j), and angular 
momentum operators Lz{j) = a_^{j)a+{j) — d_{j)a^{j), 
L±{j) — a\.{j)a^{j)/2. From the properties of the har- 
monic oscillator wave functions one finds that in the HA 
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FIG. 2: (color on-line) (a) Filling and hopping of the Px (red 
line) and Py (blue line) orbitals for u = 3/2 obtained with 
ED method on a lattice with L — 6 sites. Results agree 
with corresponding results obtained for L — 4 (thick black 
lines) and DMRG calculations (not showed since practically 
not distinguishable from ED data), (b) The energy difference 
between the two ground states in even and odd subspaces of 
the eigenstates of the symmetry operator <S. The energies are 
obtained with the ED method on the lattice with L = 4, 6, and 
8 sites (thin black, dashed blue, and solid line respectively). 
Note that corresponding lines cross the zero energy L times, 
(c) Expectation value of the staggered angular momentum 
jCz/L as a function of tunneling obtained with DMRG on 
the lattice with L = 64 sites. Nonvanishing value of C^ is 
present only in the region where ground state is degenerate. 
In all figures the shaded region denotes the range of tunnelings 
where the ground state of the system is degenerate in the 
thermodynamic limit. 



X = 6 = for any lattice configuration. In this limit 
[H{j),Lz{j)] = 0, eigenvalues of Lz{j) become a good 
quantum numbers and one recovers the situation in ;2C| . 
But for optical lattice potential (|31) we find that X,6 ^ 0, 
and consequently [H{j),Lz{j)] ^ breaking the local 
axial symmetry. 

The properties of the ground state crucially depend on 
contact interactions. As an example we revisit the case of 
total filling v = 3/2. For chosen lattice parameters, with 
proper values of contact interactions, we perform exact 
diagonalization (ED) in full many-body basis on the lat- 
tice with L — 4, 6, and 8 sites and periodic boundary 
conditions. In HA we find that for small tunnelings the 
GS of the system is degenerate, i.e. both ground states 
|Gei,en) and |Godd) have the same energy. In this way we 
reproduce the results of I20[| with all their consequences. 



However, when the anharmonicity is taken into account 
the GS looses its degeneracy and only one of these states 
becomes the ground state of the system (Fig. [1^): in 
the limit of small tunnelings GS reaches the insulating 
state in py orbital with one boson per site, and fractional 
super-fluid state in p^ orbital. Moreover, we do not find 
any significant correlations {o.l.{j)ay{j)) in this limit. In 
contrast we find that for large enough tunneling all parti- 
cles occupy Px orbital in superfluid phase (large hopping 
ha = j;^j{0'a{j)0'a{j + I)))- Thcsc results were con- 
firmed for larger L = 64 system using density matrix 
renormalization group (DMRG) approach 22, [23 1 . 

The most interesting scenario is realized for inter- 
medite values of the tunneling. As it is visible in Fig. 
[T^, the balance N^ — Ny has opposite sign in two limit- 
ing cases {tx — >■ and \tx\ — >• c«). Therefore, there ex- 
ists a particular tunneling value for which both orbitals 
are balanced. To check this point more carefully we plot 
the energy difference between ground states Eodd — Eeven 
as a function of tunneling for different lattice sizes (Fig. 
[2b)- We find that near the balanced tunneling point both 
ground states are degenerate. In fact, increasing a lattice 
size L, we find that the degeneracy occures exactly for 
L different values of the tunneling within a certain finite 
range. The range of tunneling for which Eodd — Egyen — 
does not grow with a lattice size, but saturates. With this 
observation we claim that in the thermodynamic limit 
the degeneracy of the ground state is dynamically recov- 
ered in a certain well defined range of tunnelings. At this 
region, whenever the tunneling is changed, system trans- 
fer one particle between orbitals to minimize the energy. 
Since there is no corresponding term in the Hamiltonian 
this transfer is directly related to the flip from one eigen- 
subspace of S to the other. 

In the region of recovered degeneracy both ground 
states |Gei,en) and \Godd) have the same energy. It is 
obvious that in principle any of their superposition is 
a ground state of the system. However, in the thermo- 



dynamic limit, due to the einselection principle 2J], the 
macroscopic state which is realized physically should ex- 
hibit as low entanglement as possible. To find the su- 
perposition fulfilling the lowest entanglement assump- 
tion we look for such a combination of ground states 
in which an entanglement entropy for one lattice site 
is the lowest. We minimize the von Neumann entropy 
defined as S{9, </?) = — J^i ^i l^i ^i as a function of an- 
gles 9 and (p, where A's are the eigenvalues of the re- 
duced density matrix of a single lattice site. With 
this procedure we find two orthogonal ground states 
\'^±} — (|Ge^en)±*|Godd))/'\/2 with the lowest entropy. For 
these states the reduced density matrix has two dominant 
eigenvalues Ai = A2 ^ 1/2 which only insignificantly vary 
with the tunneling. The same two states minimize the 
von Neumann entropy of the subsystem of two lattice 
sites. Then the reduced density matrix has three domi- 
nant eigenvalues Ai = A2 ^ 1/3, and A3 ~ 1/6. 



The situation is very similar to the situation in the 
standard Ising system in transverse field. In that case the 
system has also additional symmetry of flipping all spins. 
After diagonalization of the Hamiltonian in two indepen- 
dent eigen-subspaces of the symmetry operator one finds 
two degenerated ground states |±) = (|UP) ± |D0WN))/y2, 
where |UP) and JDOWN) denotes states with all spins 
up/down respectively. Obviously, from the model point 
of view any superposition of these states is a ground 
state of the system. In the thermodynamic limit, due 
to the einselection principle, one finds that the symme- 
try of the ground state is spontaneously broken and only 
|UP) or |D0WN) state can be physically realized. From this 
perspective our system has also Ising-like Z2 symmetry 
connected to the symmetry operator S. This symme- 
try is spontaneously broken in the thermodynamic limit 
to one of the two states |G±). Note, however, that in 
contrast to the usual Ising model, in the studied case 
the broken symmetry states are complex superpositions 
of original ground states. It is very unusual situation 
since the original Hamiltonian in the Cartesian basis, as 
well as in the angular momentum basis, is represented 
by a purely real matrix. Complexity appears only due to 
the additional assumption that 'Schrodinger cat' states 
can not be obtained in the thermodynamic limit. This 
assumption has many nontrivial consequences for our 
system. One of them is observation that macroscopi- 
cally realized ground state has non vanishing correlations 
Ca/sU) — (fiJ»(j)o/3(j)) for a 7^ /3 (as a consequence of 
the superposition of states with different Nx and Ny). In 
addition, the non trivial complex factor in 'superposed' 
state is responsible for the sign inversion symmetry of the 
correlations Cxy{j) — ~Cyx{j)- These leads directly to 
the observation that the broken symmetry ground state 
violates the time-reversal symmetry since this is the state 
with non vanishing expectation value of staggered angu- 
lar momentum operator tz = J2ji~^y ^z{j) (Fig- Efc)- 

To conclude, we show that anharmonicity of the op- 
tical lattice sites plays a crucial role when one studies 
the orbital properties of the system. The assumption of 
harmonicity of lattice sites is highly oversimplified since 
it leads directly to the rotational invariance of the local 
Hamiltonian in all possible arrangements of the optical 
lattice. We show that in asymmetric rectangular lattices 
it is possible to obtain degeneracy between single-particle 
energies of p-orbitals by a proper adjustment of the lat- 
tice parameters. However, this single particle levels' de- 
generacy is always lifted by anharmonicity when contact 
interactions are taken into account. Instead we iden- 
tify another symmetry operator corresponding to eigen- 
states containing even and odd number of particles in 
specific p-orbitals. We find that the degeneracy between 
them is dynamically restored due to tunneling. Addi- 
tionally we show that in the thermodynamic limit, the 
ground state breaks time-reversal symmetry along with 
the U{1) X C/(l) symmetry down to Ising-like Z2 sym- 



metry. Examples of such dynamical restoration and then 
breaking of symmetry is quite rare in condensed mat- 
ter physics though such cases are known in high-energy 
physics (i.e. Schwinger mechanism for dynamical gener- 
ation of mass). The studied system can be probed exper- 
imentally by trapping ultracold atoms in optical lattice 
setup. 
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